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We present a formulation of collisional gyrokinetic theory with exact conservation laws for energy 
and canonical toroidal momentum. Collisions are accounted for by a nonlinear gyrokinetic Landau 
operator. Gyroaveraging and linearization do not destroy the operator’s conservation properties. 
Just as in ordinary kinetic theory, the conservation laws for collisional gyrokinetic theory are selected 
by the limiting collisionless gyrokinetic theory. 


Introduction. — One of the greatest unsolved problems 
in the theory of magnetically-confined plasmas is under¬ 
standing and controlling the turbulent flux of particles 
and heat into a fusion reactor’s wall [1]. It is believed that 
the predominant cause of these fluxes is low-frequency 
fluctuating electromagnetic fields with wavelengths on 
the order of the gyroradius. While a collisionless gy¬ 
rokinetic model of these fluctuating fields has been de¬ 
veloped that is fully consistent with the First Law of 
Thermodynamics (for a recent review see Ref. [2]), this 
energetically-consistent model has the serious flaw of ig¬ 
noring collisions altogether. 

In order to accurately describe irreversible plasma 
transport processes, the effects of collisions must be in¬ 
corporated into gyrokinetic theory. Previous work on lin¬ 
ear gyrokinetic collision operators [3-5] assumed a strict 
two-scale separation between a large-scale equilibrium 
distribution function Fo and a small-scale fluctuating 
part SF = F — Fo. Conservation properties of the colli¬ 
sion operator in Ref. [3], for example, were discussed in 
the gyroBohm limit. Here, we will focus on nonlinear gy¬ 
rokinetic collision operators for a global full-F approach 
that do not make this split, and that can thus investigate 
more completely the possible effects of finite e = Pi/L 
in experiments, such as corrections to gyroBohm scaling 
and non-local turbulence spreading (see footnote 5 on p. 
427 in Ref. [2].) 

When flnite-e effects are accounted for, preserving ex¬ 
act conservation properties, and therefore ensuring con¬ 
sistency with the First Law of Thermodynamics, is a 
nontrivial unsolved problem. The collision operators in 
Refs. [3, 4], for example, were obtained by transforming 
a particle-space collision operator with exact conserva¬ 
tion properties into the lowest-order guiding center co¬ 
ordinates. While this approach guarantees the existence 
of energy and momentum-like quantities that annihilate 
the collision operator, these same quantities are not con¬ 
served by the full-F collisionless gyrokinetic system, and 


therefore fail to be conserved by the full-F collisional sys¬ 
tem. More generally, existing gyrokinetic collision opera¬ 
tors are not energetically consistent in a full-F formalism 
because: (a) the gyrocenter coordinate transformation, 
and therefore any collision operator transformed into gy¬ 
rocenter coordinates, is only known as an asymptotic ex¬ 
pansion in the gyrokinetic ordering parameter e; and (b) 
replacing the asymptotic expansion of such an operator 
with a truncated power series destroys exact conserva¬ 
tion laws. The purpose of this Letter is to present the 
first collisional formulation of global full-F gyrokinetics 
with exact conservation laws. 

Electrostatic Model. — For the sake of simplicity, our 
discussion will focus on quasi-neutral electrostatic gyroki¬ 
netics (for instance, see Ref. [6]). However, the ideas be¬ 
hind our discussion apply equally-well to electromagnetic 
gyrokinetics (for example, see Ref. [7].) Our primary re¬ 
sult consists of an expression for the non-linear Landau 
operator in gyrocenter coordinates that is corrected by 
small terms to ensure exact energy and momentum con¬ 
servation [see Eq. (24).] These correction terms are analo¬ 
gous to the Hji-denominators in the Hamiltonian guiding 
center theory introduced by Littlejohn [8]; they do not 
increase the theory’s order of accuracy, but they are es¬ 
sential to include for the sake of ensuring exact energy 
and momentum conservation. 

As a first step, we review how the energy conserva¬ 
tion law is discussed in collisionless kinetic theory. The 
governing equations of collisionless electrostatic kinetic 
theory are the Vlasov-Poisson equations, 

dtfs + {fs,Hs}=0 ( 1 ) 

= -47rp(/), (2) 

where fa is the species-s distribution function, ip is the 
electrostatic potential, p{f) is the charge density, Hg = 
p'^ f2ms+esP, and {•, •} is the standard canonical Poisson 
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bracket. Equations (l)-(2) conserve the total energy 

£ = + + p) 

where (•, •) denotes the standard E^-pairing of functions 
on configuration space and dz = dx dp. Because binary 
collisions conserve energy, Eq. (3) must also be conserved 
in collisional kinetic theory. In particular, if the Vlasov- 
Poisson equations are modified by the addition of a bi¬ 
linear collision operator, 

dtfs + {/., H,} = ^ Cs-s{fs, fs) ( 4 ) 

S 

A(p =-Anpif), (5) 

then Css must be chosen to satisfy the condition 

® = 

= / HsCssifs,fs)dz. ( 6 ) 

Because this identity must hold for an arbitrary multi¬ 
species distribution function, the collision operator there¬ 
fore has to satisfy the well-known identities 

J HsCssifsJs) dz + J HsCssifs, fs) dz = 0 , (7) 

which express the fact that the energy gained by species 
s due to collisions with species s is precisely the energy 
lost by species s due to collisions with species s. The 
non-linear Landau operator (summation rule is implied), 

Css{fsJs) = -^{Xi,jn, ( 8 ) 

satisfies the identities (7), and therefore defines an 
energetically-consistent collisional kinetic theory. Here 
Tss = 47rege?lnA; the 3-component vector 7®'* is 

lt{z) = j 5{x- x) Q®®(z, z) Ass{z, z) dz\ (9) 

the 3x3 matrix Q®® is given by 

L (10) 

Wssyz^z) 

where P(^) = I — is the orthogonal projection onto 
the plane perpendicular to the vector the velocity dif¬ 
ference Wss is given by 

Wss{z,z) = {x,Hs}(z) - {x,Hs}{zy, (11) 

and the vector 


When comparing this form of the Landau operator to 
more conventional expressions, it is useful to note that 
{x,g} = dpQ, where g is any function on phase space, 
so that the collision operator (8) describes collisions in 
momentum space. Moreover, the identities (7) follow im¬ 
mediately from the fact that the velocity difference Wss 
is a null-eigenvector of the matrix Q®®. 

Electrostatic Gyrokinetic Model. — In order to apply 
this same argument to gyrokinetic theory, we start with 
the gyrokinetic Vlasov-Poisson system 

dtFs + {Fs,Hiy}f = 0 (13) 

y-p = p{F). (14) 

Here, Fs is the gyrocenter distribution function; ip is the 
electrostatic potential; {•, •}f'^ is the guiding center Pois¬ 
son bracket; 

Hfy = iJf + e. W + |({^, d-lf) ^ Ks{E) + csp 

(15) 

is the gyrocenter Hamiltonian; ip{z) = p{X -|- Pos), 
where Pos is the lowest-order guiding-center gyroradius; 
{•) denotes the gyroaverage; if denotes the gyroangle an¬ 
tiderivative oi Ip = Ip — (ip); Ks{E) is the gyrocenter 
kinetic energy; P = —6IC/6E is the gyrocenter polar¬ 
ization density; 1C = Jf,sl FsKs{E) dz^^\ and dzf^ de¬ 
notes the guiding center Liouville volume element. These 
equations govern collisionless quasineutral electrostatic 
gyrokinetic theory in the “high-flow” regime (see [9] and 
references therein) and they conserve the total energy, 

£^y = FsHfydzr, (16) 

s 

exactly. Note that the quasineutrality equation (14) im¬ 
plies that this system governs plasma dynamics on time 
scales long compared to the period of plasma oscillations. 

The equations governing collisional gyrokinetic theory 
are given by adding a bilinear collision operator to the 
gyrokinetic Vlasov-Poisson equations, 

dtFs + {Fs, Elf }f = ^ Cfl {Fs,Fs) (17) 

S 

V-P = p{F). (18) 

Because the conservation laws of ordinary collisional ki¬ 
netic theory are consistent with those of collisionless ki¬ 
netic theory, the gyrokinetic collision operator Cff must 
not alter the conservation of £^. Thus, 

0 = E / H^dtPs dzf + (p(F) - V • P, dtp) 

= J2[ HiycfUFs,Fs)dzic 


Ass{z,z) = fs{z){x,fs}{z) - {x,fs}{z)fs{z) 


( 12 ) 


(19) 
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This identity will be satisfied for a general multi-species 
gyrocenter distribution function if and only if 

I HfCf, (F,, F,) dzf + I Hfdl{F,, F,) dzf = 0, 

( 20 ) 

which is the gyrokinetic version of Eq. (7). The identi¬ 
ties (20) must be satisfied exactly by any energetically- 
consistent gyrokinetic collision operator. 

An energetically-consistent collision operator — While 
Eq. (20) imposes important qualitative constraints, they 
cannot determine the form of the gyrokinetic collision op¬ 
erator by themselves. A quantitative constraint is nec¬ 
essary as well. To this end, it is important that the gy¬ 
rokinetic collision operator agrees with the the transfor¬ 
mation of the particle-space Landau operator [10] into 
gyrocenter coordinates, at least up to some desired order 
in the gyrokinetic ordering parameter e. Is it possible 
to satisfy these qualitative and quantitative constraints 
simultaneously? The answer is “yes”. 

We have discovered an accurate gyrokinetic collision 
operator that is consistent with the conservation laws 
of collisionless gyrokinetic theory, and therefore the first 
law of thermodynamics. The form of the operator is 
suggested by the somewhat-peculiar presentation of the 
particle-space Landau operator given earlier. Let j/s = 
X + pos and define the gyrocenter velocity difference 

W!i{z, z) = {y«, Ff }f (z) - {y„ Ff jf (z), (21) 

the associated 3x3 matrix 

^-)], (22) 

and the vector 

Afs{z,z) = Fs{z){ys,Fs}f{z) - {Vs, Fs}T{z)Fs{z). 

(23) 

The energetically-consistent gyrokinetic Landau operator 
is given by 

dUFs,F,) = , (24) 

where 

'l'gy(^) = J z)Qlyiz, z)Afg{z, z) dzf, (25) 

and Sf{z,z) = 5{ys{z) — ys{z)). Note that this oper¬ 
ator depends explicitly on the electric field through the 
gyrocenter Hamiltonians that appear in Eq. (21). Using 
a straightforward, but tedious argument that is not re¬ 
produced here, we have shown that this operator agrees 
with the Landau operator transformed into gyrocenter 
coordinates with leading-order accuracy. 


Because the proof is simple, we will now show explicitly 
that the gyrokinetic Landau-Poisson system (17) defined 
in terms of the collision operator (24) has exact conser¬ 
vation laws for energy and momentum. We hope to con¬ 
vey the similarity of this demonstration with the analo¬ 
gous demonstration for the ordinary Landau-Poisson sys¬ 
tem (4)-(5). However, a word of caution is in order here. 
It is essential that the guiding center Poisson brackets 
that appear in Eq. (24) be genuine Poisson brackets (i.e., 
the brackets must satisfy the Leibniz and Jacobi identi¬ 
ties). Dropping terms from a bracket that satisfies these 
properties will destroy the gyrokinetic Landau-Poisson 
system’s exact conservation laws. 

Energy conservation — Proving that the gyrokinetic 
Landau operator (24) satisfies the identities (20) is very 
similar to proving that the particle-space Landau op¬ 
erator satisfies the identities (7). Setting Egg = 
f F[f'Cgs{Fg,Fs) dzf, it is simple to verify that 

fs. + ^-ss = ^ JJ (WS')tQI® dzf dzf, (26) 

where all two-point quantities in the integrand are eval¬ 
uated at {z,z) and 4 denotes the ordinary matrix trans¬ 
pose. Because Q®® is a symmetric matrix with null eigen¬ 
vector W®/, the right-hand-side of this equation van¬ 
ishes exactly. Thus the gyrokinetic Landau operator (24) 
satisfies the identities (20) exactly, and the gyrokinetic 
Landau-Poisson system (17) has an exact energy conser¬ 
vation law, d£^^/dt = 0. 

Toroidal momentum conservation — We will prove 
that if the background magnetic field is axisymmetric, 
then the gyrokinetic Landau-Poisson system conserves 
the total toroidal momentum 

= E (27) 

s -7 

where p^s is the guiding center canonical toroidal mo¬ 
mentum [11]. If the background magnetic field has ad¬ 
ditional symmetries, a similar proof of the conservation 
of the corresponding total momentum can easily be con¬ 
structed. The time derivative of Eq. (27) yields 

^ = E / Pc^sC!UFg,F-g) dzf = hss, (28) 

s,s s,s 

where is conserved exactly by the gyrokinetic Vlasov- 
Poisson system. Here, we find 

P<pSS T P(pSS — 

(29) 

Now using the fact that p^s is the generator of infinites¬ 
imal toroidal rotations, we can see that {ys,P(j,s}f = 
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e-z X Vs, where is the unit vector along the axis of 
rotation. Therefore the vector quantity {{ysTP4>s}T ~ 
{ys,P<j,s}f)S^^ = X {ys - ys)S^1 = 0, which follows 
from standard (5-function properties. This shows that 
P(j>ss + Pcpss = 0, which in turn implies total toroidal 
momentum conservation dPrf,/dt = 0. 

Entropy production — As we have discussed, these 
conservation laws ensure that the gyrokinetic Landau- 
Poisson system is consistent with the the First Law of 
Thermodynamics. On the other hand, they do not di¬ 
rectly imply that the gyrokinetic Landau-Poisson system 
is consistent with the Second Law of Thermodynamics. 
To verify that entropy is indeed a non-decreasing func¬ 
tion of time, we have computed the time derivative of 
S = Ps^nFg dzf‘^ and found 

f = ^ // d-zf dzf. (30) 

Because Qg® is a positive semi-definite matrix and the 
distribution function is positive [12], the right-side of 
Eq. (30) is non-negative, which is the desired result. 

Note that this proves one “half’ of a gyrokinetic ver¬ 
sion of Boltzmann’s iL-theorem. The missing ingredient 
is a complete characterization of the distributions that 
satisfy dS/dt = 0, i.e. the gyrokinetic Maxwellians. Be¬ 
cause the guiding center Poisson bracket is rather compli¬ 
cated, we have not yet found a complete characterization. 
However, we have verified that the distribution 

1 / Rsy \ 

Fms = — exp ^ , (31) 

where Zg = J exp(—iLf^/T) dzf^ is the partition func¬ 
tion, maximizes the entropy. We leave the characteriza¬ 
tion of the most general gyrokinetic Maxwellian, which 
would be useful for the sake of deriving dissipative gy- 
rofluid models with exact conservation laws [13], as a 
topic for future study. 

Gyroaveraging — When the collision frequency is much 
smaller than the gyrofrequency [14], the full gyrokinetic 
Landau operator (24) can be replaced with that opera¬ 
tor’s gyroaverage, (Cff). When this is done, the gyroki¬ 
netic Landau-Poisson system becomes the gyroaveraged 
Landau-Poisson system, 

dtPs + {Fs, {Fs,F-g)) (32) 

S 

V • P = P{F), (33) 

where Fg is now interpreted as the gyroaveraged part of 
the distribution function. Because the functions and 
p^g are independent of the gyrophase, the proofs of en¬ 
ergy and momentum conservation given earlier work with 
Cgg replaced by (Cff). Thus, the gyroaveraged Landau- 
Poisson system has exact energy and momentum conser¬ 
vation laws. 


Linearization — Closely related to the gyroaveraged 
Landau-Poisson system is the collisionally-linear gyroav¬ 
eraged Landau-Poisson system, 

dtFg + {Fg, Hfy}f = ^ (^6Cir + , (34) 

s 

V • P = p{F), (35) 

where the linearized test-particle and field-particle colli¬ 
sion operators are 

SCiriFg) = (d (P., Pm.-)), (36) 

SCff'^iFg) = (Cf|(PM«,P.)). (37) 

This system of equations is obtained from the gyroaver¬ 
aged Landau-Poisson system by assuming Fg = Fms + 
SFg and then dropping the non-linear term in the collision 
operator, {SFg, SFg)). Note that (Fms, Fms)) = 
0 [15]. Because the gyrokinetic Landau operator sat¬ 
isfies the identities (20), it is straightforward to prove 
that these equations have the same conservation laws 
for energy and momentum as the gyroaveraged Landau- 
Poisson system. 

Concluding remarks — The key to deriving an 
energetically-consistent formulation of collisional gyroki- 
netics was first expressing the particle-space Landau op¬ 
erator in terms of Poisson brackets “as much as possi¬ 
ble,” which was an idea first championed by Brizard in 
Ref. [14]. In particular, the identity 

V - V = {x,Hg}{z) - {x,Hg}{z) (38) 

suggests that the appropriate definition of the gyrocen- 
ter velocity difference is given by Eq. (21). This idea, 
together with the procedure given earlier for determin¬ 
ing the energetic consistency constraints, appears to be 
appropriate for deriving energetically-consistent collision 
operators for other reduced plasma models as well. In 
future work, we will report on the energy-conserving col¬ 
lisional formulations of electromagnetic gyrokinetics and 
oscillation center theory. 

We note that, although the gyrokinetic Landau opera¬ 
tor (24) and its linearized forms (36)-(37) may prove dif¬ 
ficult to implement numerically, they identify the proper 
formalism for the inclusion of collisional transport in gy¬ 
rokinetic theory. Hence, these gyrokinetic collision oper¬ 
ators form the basis from which approximations can be 
implemented for practical applications. 

Lastly, by setting (/? = 0 in the above formulas, our re¬ 
sults reduce to an energy-momentum-conserving guiding 
center collision operator. This operator would be ide¬ 
ally suited to incorporating collisions into orbit-following 
codes such as ORBIT [16]; see Ref. [17] for recent work 
on the Monte Carlo implementation of a 5D guiding cen¬ 
ter Fokker-Planck collision operator. All previous guid¬ 
ing center collision operators that have been applied in 
orbit-following codes either resort to ad hoc methods to 
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ensure exact conservation laws [18], or else inconsistently 
account for inhomogeneities in the magnetic field [19]. 
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